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Abstract. We present the far-UV spectrum of the quasar HS 1700+6416 taken with FUSE. This QSO provides the second 
line of sight with the Hen absorption resolved into a hya forest structure. Since HS 1700+6416 is slightly less redshifted 
(z em = 2.72) than HE 2347-4342, we only probe the post-reionization phase of He n, seen in the evolution of the He n opacity, 
which is consistent with a simple power law. The Hen/Hi ratio 77 is estimated using a line profile-fitting procedure and an 
apparent optical depth approach, respectively. The expected metal line absorption in the far-UV is taken into account as well as 
molecular absorption of galactic H 2 . About 27 % of the 77 values are affected by metal line absorption. In order to investigate the 
applicability of the analysis methods, we create simple artificial spectra based on the statistical properties of the H 1 hya forest. 
The analysis of the artificial data demonstrates that the apparent optical depth method as well as the line profile-fitting procedure 
lead to confident results for restricted data samples only (0.01 < T Hr < 0.1 and 12.0 < logiVni < 13.0, respectively). The reasons 
are saturation in the case of the apparent optical depth and thermal line widths in the case of the profile fits. Furthermore, 
applying the methods to the unrestricted data set may mimic a correlation between the Hen/Hi ratio and the strength of the 
Hi absorption. For the restricted data samples a scatter of 10 — 15% in 77 would be expected even if the underlying value 
is constant. The observed scatter is significantly larger than expected, indicating that the intergalactic radiation background 
is indeed fluctuating. In the redshift range 2.58 < z < 2.72, where the data quality is best, we find n ~ 100, suggesting a 
contribution of soft sources like galaxies to the UV background. 
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1. Introduction 

The shape and the strength of the intergalactic UV background 
plays an important role in governing the evolution of the ion- 
ization state of the intergalactic medium (IGM). Since the IGM 
is observed to be highly ionized, ionization corrections have 
to be applied to measure properties of interest like its metal- 
licity. The strength of the UV background, at least at the Hi 
ioniz ation edge, can be derived from the proximity effect (e.g. 
IScott et afll2000l) . whereas more information is needed to con- 
Send offprint requests to: C. Fechner 



strain the shape of the ionizing radiation field. Agaf onova et alJ 
(2005) introduced a method to estimate the shape of the ioniz- 
ing background from the observation of metal line absorbers. 

Theoretical calculations model the shape of the UV back- 
ground as the radiation of sources attenuated while propagat- 
ing though the IGM. Quasars have been adopted as sources 
jHaardt & Madaull99fjlFardal et all 1991 whose radiation has 
been filtered by the IGM. Since hydrogen and helium are the 
most abundant elements in the IGM, the resulting shape of 
the mean UV background is characterized by a sharp break at 
1 Ryd and another break at 4 Ryd. 
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Observational results, however, suggest that the intergalac- 
tic UV background is dominated not only by the filtered ra- 
diation of quasars but also by a significant contribution from 
the softer radiation of galaxies. This evidence come s from 
the study of the Hi opacity ("e.g. iKirkman et aDl2005h or the 
number density evo l ution of Hi absorbers (e.g. Zhang etafl 
1 19971 iBianchiet alJ 1200 ll) . and also the derivation of rea- 
sonable metallicities from observed metal line absorption re- 

i II 1 

quires softer radiation ( Aguirre et al. 2004). Recent computa- 
tions of the UV background therefore consider the contribution 
of galaxy radi ation ("e.g.lHa ardt & Madaul200ll) . Furthermore, 
ISchavd (12004) and Miralda-Escude (2005) pointed out that lo- 
cal sources may dominate over the mean background. 

Besides Hi opacity measurements and the investigation 
of metal line systems, Hen provides another observational 
probe of the UV background. Comparing the strength of the 
He ii Lya absorption to that of Hi leads to the column den- 
sity ratio rj = Nnen/NHi, which is directly related to the ion- 
izing backgroun d at th e Hi and Hen photoionization edge 
(e.g. Far dal et alJll998i) . Low values are expected for hard ra- 
diation, while high values should be measured if the ioniz- 
ing radiation is dominated by soft sources. However, the ob- 
servation of the He ii Lya transition (A lest = 303.7822 A) is 
very difficult. An unabsorbed, UV-bright QSO is needed at 
redshift z > 2 to be ob servable with HST/STIS or FUSE 
JPicard & .Takobse^ 19931). Up to now six of such quasars 



thors (e.g .[Reimers et al. 1992: Vogel & Reimers II 19931 


1995; 


Kohleret aO fl996: 


Petitiean et al. 1996; Tripp et al. 


1997; 


Simcoe etall2002. 


2006). Thus, metal line absorption features 
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have been found ( Jakobsen et al. 1994; Davidsen et al. 1996; 
Reimers et al. 1997; Anders onet al.lll999HZheng et ai]E( )04a; 



Rei mers et al] l2005al) . Comparisons of the Hen opacity to 
the Hi opacity along a number of these sightlines point to- 
ward a dominant contribution to the ionizing UV background 
from hard sources, such as quasars, but with significant con- 
tributions also from s o ft sources such as star-forming galax- 
ies .Heap etaD 12000. iKriss et all 1200 it ISmette et all .2002; 
IZheng et al.l2004aUShull et al 120041 

The He ii Lye forest has been resolved fo r the first time to- 
wards the quasar HE 2347-4342 with FUSE dKriss etall2001_ 
IShuU et al.l2004HZheng et al.l2004bl) . Different analysis meth- 
ods have been applied to derive the column d ensi ty ratio (line 
profile f itting was used bvlKriss et alJ d200ll) and lZheng et alJ 
(l2004bll . while IShull et al.1 _2004l) applied an apparent optical 
depth method) and lead to the same results: The UV back- 
ground is highly variable on small scales with a fluctuation 
of the column dens ity ratio be t ween 1 < 77 < 1000 or even 
more. Additionally. IShull et al.1 -2004) find an apparent corre- 
lation between the 77 value and the strength of the H 1 absorption 
in the sense, that 77 appears to be larger in voids. 

In this work we present FUSE observations of the quasar 
HS 1700+6416, which provides the second line of sight of 
a resolved Hen Lya forest. Owing to the v ariability of the 
UV flux of this QSO .Reimers et alJ l2005bl) . the data is of 
comparable quality like those of HE 2347-4342 (S /N ~ 5, 
R ~ 20000). The emission line redshift of HS 1700+6416 is 
Zem = 2.72, i.e. this line of sight probes the post-reionization 
universe, while in the spectrum of HE 2347-4 342 the end of 
the H e 11 reionization phase has been detected .Reimers et alJ 
Il997l) . HS 1700+6416 has a rich metal line spectrum in 
the optical as well as in the UV comprising seven Lyman 



are also expected to arise in the far-UV polluting the Lya- for- 
est absorption. Therefore, we will include a model of the metal 
line absorption in the analysis to obtain unbiased results. For 
this purpose, we derived a pre diction of the metal line spectrum 
IFechner et al. 2006), which we will consider in 



in the far-UV 1 
this work. 

This paper is organized as follows: After presenting the ob- 
servations in Sect. [5] and describing the extrapolation of the 
continuum for the FUSE spectral range in Sect. |3] we present 
the fa r-UV metal line prediction derived in IFechner et alJ 
(2006) in Sect. 0] and give a comparison to the data. In Sect. 
|5]we discuss our artificial data and analyze them using the line 
profile-fitting procedure and the apparent optical depth method 
to investigate the applicability of these methods. The results 
derived from the observed data are presented and discussed in 
Sect. [6] We conclude in Sect.0 

2. Observations 

The Far Ultraviolet Spectroscopic Explorer (FUSE) uses four 
independent optical channels to enable high-resolution spec- 
troscopy in the far-ultraviolet wavelength range below 1200 A. 
In each channel a primary mirror gathers light for a Rowland- 
circle spectrograph. Two-dimensional photon-counting detec- 
tors record the dispersed spectra. Two of the optical channels 
use LiF coatings on the optics to cover the 990 - 1 187 A wave- 
length range. The other two channels cover shorter wavelen gths 
down to 912 A using SiC-coated optics. |Moos et a"i1_20 00) give 
a full description of FUSE, and ISahnowetaiT J20001) summa- 
rize its in-flight performance. 

We observed HS 1700+6416 with FUSE through the 30"- 
square low-resolution apertures (LWRS) during four epochs 
in October 1999, May 2002, February/March 2003, and May 
2003. Table ^ gives details of the individual observations. 
Between the exposures in May 2002 and 2003 the UV bright- 
ness of HS 1700+6416 increased by a factor ~ 3 (see 
iReimers et al.ll2005bll . Therefore, the final data are of better 
quality than we originally expected. The total exposure time 
is roughly 791 118 s with 290 230 s during orbital night. 

We use the FUSE calibration pipeline CALFUSE V2.2.1 
to process the data. The earlies t version s of th e calibration 
pipeline are described by Sahno w et al ] i200(j . CALFUSE 
V2.2.1 is described in detail on the FUSE web site 
(http://fuse.pha.jhu.edu/analysis/pipeline_reference.html). The 
major improvements to the pipeline since the description of 
Sahn ow et alJ .20001) that are relevant for our work are that 
data screening for bursts has been automated, and that a two- 
dimensional, multicomponent model for the background has 
been implemented. The background model takes into account 
the uniform (but time-varying) particle background in the de- 
tector as well as the spatially varying and time-varying portions 
of the background due to scattered light within the instrument. 

Since we are working at the very limits of instrument sen- 
sitivity, we take additional non-standard steps in our process- 
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Table 1. Log of FUSE Observations of HS 1700+64 



Observation ID 


Date 




'total (s) 


'night (S) 


PI 100401 


1999 


Oct 


14 


57 


038 


23 275 


C1230101 


2002 


May 


15 


54 


710 


17 702 


C1230102 


2002 


May 


16 


53 


821 


15 663 


C1230103 


2003 


Feb 


27 


156 


963 


93 828 


C1230104 


2003 


Mar 


03 


130 


605 


51 584 


P3060101 


2003 


May 


02 


61 


326 


18 520 


P3060102 


2003 


May 


03 


57 


952 


17 863 


P3060103 


2003 


May 


10 


91 


590 


23 824 


Z0120101 


2003 


May 


03 


53 


896 


14943 


Z0120102 


2003 


May 


04 


43 


520 


9702 


Z0120103 


2003 


May 


05 


21 


064 


3019 


Z0120104 


2003 


May 


10 


9 


916 


307 



ing to reduce the background. We restrict the data to orbital 
night only (scattered light is ten times brighter during orbital 
day), and we select only detector events with pulse heights 
in the range 4 - 16. This reduces the background by an addi- 
tional factor of two while still retaining photometric accuracy 
in the final result. Finally, since the instrument throughput is 
roughly a factor of three higher in the LiF channels compared 
to the SiC channels, we only use data from the LiF channels. 
The final spectrum therefore only covers the wavelength range 
1000 - 1 180 A, with a gap at 1080 - 1086 A due to the detector 
layout. 

Each observation interval listed in Tablets processed sep- 
arately. To assure the best spectral resolution and an accurate 
wavelength scale, we cross-correlate the overlapping wave- 
length regions of each extracted spectrum with each other. No 
wavelength shifts were identified that exceeded our final bin 
size of 0.05A. The extracted spectra from each detector seg- 
ment and exposure are then merged onto a uniform wavelength 
scale in 0.05A bins. To maximize the resulting signal-to-noise 
ratio, each extracted spectrum is weighted by the product of 
the exposure time and the continuum flux at 1000 A. This ef- 
fectively weights each spectrum by the expected relative num- 
ber of photons in each spectrum, which is the optimal weight- 
ing for data with Poisson-distributed errors. The resulting spec- 
trum shown in Fig. ^ has S /N ~ 5 and a spectral resolution of 
20000. 

In the final merged spectrum we notice some potential 
problems with the zero flux level, particularly in the data at 
wavelengths shortward of 1080 A. Strong, saturated H i absorp- 
tion is detected at z ~ 2.315 and z ~ 2.433 arising from LLS 
at these redshifts. At the corresponding positions in the FUSE 
spectrum (~ 1007 A and ~ 1043 A) strong but slightly unsat- 
urated features are detected. Since the column density ratio 77 
should be > 1 because of physical plausibility, these features 
actually have to be saturated. However, the troughs of these 
presumably saturated lines are consistent with zero flux within 
the errors of our S /N ~ 5 spectrum. We note the systematic er- 
rors that may result from this in our future discussion, but note 
that data at longer wavelengths (A > 1086 A) does not appear 
to have the same problem. Saturated Hi features at z ~ 2.597, 



2.602, 2.634, and 2.683 show also saturated Hen counterparts 
at ~ 1092 A, 1094 A, 1104 A, and 1119 A in the FUSE spec- 
trum (see e.g. Fig [3}. Thus, we are confident we can obtain 
unbiased results in the range 2.58 < z < 2.72. Problems with 
the interpretation of the results due to inaccuracies in the zero 
flux level in the range 2.29 < z < 2.56 will be discussed later. 

The corresponding H 1 Lya forest was observed w ith Keck. 
We ha ve ac cess to two da t asets a lready published by Songaila 
dl998h and ISimcoe et alJ l2QQ2\i . These two datasets are co- 
added obtaining a resulting spectrum with a total exposure time 
of 84 200 s and a signal-to-noise of S /N ~ 100 at 4000 A. The 
co-added spectrum covers the wavelength range 3680 - 5880 A 
with a resolution of R ~ 38 500. The Simcoe data covers 
wavelengths down to ~ 3220 A with less signal-to-noise. We 
also have used the absorption lines identified in this portion of 
the spectrum to co nstrain the metal line system models (see 
iFechner et all2006l) . 

3. Continuum definition 

In order to estimate the continuum in the FUSE portion of the 
spectrum we use low resolution data taken with HST/STIS in 
May 2003 simultaneously with the FUSE observations using 
the gratings G140L and G230L. The spectra cover the wave- 
length range 1150 - 3200 A. The continuum is estimated con- 
sidering galactic extinction, Lyman limit absorption, and the 
intrinsic spectral energy distribution of the QSO. Shortwards 
of the Lyman limit break the contribution of the optical depth 
of the LLS decreases due to the v -3 dependence of the hy- 
drogen photoionization cross section. The strength of the 
break is proportion al to the H 1 column density of the systems 
(see iMol ler ffi: Jakpbser] fl99oT) . which have been measured by 
IVogel&Reimersldl995l) . 

The galactic extinction towards HS 17 00+6416 is about 
E(B - V) = 0.025 following lschlegel et alJ lll998h . whose val- 
ues are based on dust maps created from COBE and IRAS 
infr ared data. X- r ay ob servations of this quasar carried out 
by Reim erset alJ d 1995b yield a hydrogen c olumn density of 
Af HT = (2.9 + 1.3) • 10 20 cnr 2 . Following IBiplas & Savagel 

( 1994) , this corresponds to E(B - V) = 0.06. Though the latter 
value is derived from direct observations of this line of sight, 
it is based on the measureme nt of H 1 only a ssuming a galaxy- 
wide gas-to-dust ratio, while ISchlegel et all Jl998l) use a com- 
bination of H 1 and the dust distribution. Theref ore, w e adopt 
the value E(B -V) = 0.025 from Schlegel et alJ Jl998l) . 

The corrected flux is fitted with a power law f v cc v a repre- 
senting the intrinsic QSO continuum. The best fit succeeds us- 
ing two power laws with a break at 2000 A. The best fit yields 
a = -1.23 at A < 2000 A and a _ -1 54 at A > 2000 A 
using the extinction curve of lCardelli et alJ dl989h. Comparing 
the STIS spectra to the FOS data analyzed by Vogel & Reimers 

( 1995) taken in December 1991 the flux level is the same in the 
two datasets down to about 2000 A. For this wavelength range 
the spectral index is in agreement with that derived from the 
older data. The continuum rises more steeply with decreasing 
wavelength in the FOS spectrum. At 1200 A the flux level of the 
STIS data is depressed by a factor of about two in comparison 
to the FOS data. Between the two exposures in the course of 
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1000 1035 1050 1075 1100 1135 1150 1175 1300 



x(A) 

Fig. 1. The FUSE data and the extrapolated cont inuum. The solid lin e presents the continuum considering the ICardelli et alJ 
( 1989) extinction curve and E(B - V) = 0.025 (Schleg el et al.ll 1 9981). For co mpariso n we show the extr apolated continuum 
usin g E(B - V) = 0.060 adopted from X-ray observations dReimers et al! 19951) and the lCardelli et al I d 19891 dotted line) or the 
IPeil ll 19921 dashed line) extinction curve, respectively. Furthermore, we overlay the low resolution STIS data which was taken 
simultaneously with the FUSE data in May 2003. For a higher resolved presentation of the normalized FUSE spectrum see Fig. 

m 



the observations with FUSE separated about a year, the quasar 
flux increased about a factor ~ 3 in the far-UV. This means 
HS 1700+6416 is highly va riable in the intrinsic EUV even on 
relatively short time scales dReimers et all2005bl) . Thus, in or- 
der to extrapolate a reliable continuum it is very important to 
have simultaneous STIS and FUSE observations. Fig.[2illus- 
trates that the FUSE and STIS portion of the spectrum taken in 
2003 are well matched. 

A crucial point in the continuum extrapolation is the choice 
of the extinction curve especially beca use of its steep rise i n the 
UV. Adopting the analytic formula of Card elli et all dl989l) we 
find a flux level at 1000 A of ~ 2.69 • lO" 15 ergs" 1 cm 2 A 1 . 
Tests with the reddening value derived from the X-ray ob- 
servations dReimers et alJll995l) lead to a flux at 1000 A de- 
pressed by ~ 25 % (s ee Fig. [Q. Adopting the analytical ex- 
pression of IPeld 19921) an the same reddening value leads to a 
flux level slightly above the original extrapolation. But since 
the FUSE spectrum is extremely noisy especially in the lower 
wavelength range, the dominating uncertainties result from the 
profile fitting. The choice of the flux level affects the results 
only marginally. The apparent optical depth method is more 
sensitive to the continuum level. Adopting the normalization 
according to the lower continuum level given in Fig.^would 
lower the 77 values of individual bins by roughly 0.2 dex near 
1000 A. Howeve r, we are confident th at the continuum extrap- 
olation using the Carde lliet al.1 JT989) extinction curve and the 
E(B - V) from Schlege l et al.l (11998 ) represents a reasonable 
approximation of the real continuum. 

The continuum in the Keck portion of the spectrum is es- 
timated in the course of the line profile-fitting procedure. The 
line-fitting program CANDALF developed by R. Baade per- 
forms the Doppler profile fit and the continuum normalization 
simultaneously. Thus, the continuum determination in the Lyo- 
forest is more reliable in comparison to an a priori continuum 
definition. The normalized Keck spectrum in the wavelength 
range corresponding to the He n data is shown in Fig. 13 



4. Metal line absorption 

The spectrum of HS 1700+6416 is characterized by a large 
number o f metal line absorption features in the optical spec- 
tral ra nge dPetitiean et al Jl996tlTripp et a lJl997tlSimcoe et all 
| 2002l 120061) as well asliTfhe UV dReimers et al.l 1 19921: 
IVogel & Reimersll993iri995l:lKohler et alJll996l) . Since a sig- 
nificant contribution of absorption by metals is supposed to 
be present even in the far-UV, the expected metal lines in the 
FUSE portion have to b e taken into accoun t. Therefore, in the 
first paper of this series dFechner et all2006l) we construct pho- 
toionization models of the metal line systems with the main 
purpose to derive a prediction of the metal line features ex- 
pected to arise i n the FUSE sp e ctral r ange. In the following the 
main results of iFechner et al.1 d2006l) are summarized and we 
compare the predicted far-UV metal lines to the FUSE data. For 
details of the modelling procedure we refer to IFechner et alJ 
(2006). 

On the basis of the modelled column densities a metal line 
spectrum for the wavelength range 1000 - 1 180 A covered by 
FUSE is predicted. These lines (cf. Fig. 7 of IFechner et alJ 
2006) are overlayed to the observed FUSE data in Fig. [3] It 
can be seen, that some rather strong features can be identified 
as metal lines. Obviously, the prediction is consistent with the 
data for most of the lines. 

In case of the LLS at z = 0.8643, however, strong features 
of oxygen are predicted, but the observed spectrum shows less 
absorption. In particular, O iv is overestimated, as can be seen 
at ~ 1032 A a nd 1 134 A in Fig. |3] The model of this system 
is based on the lHaardt & MadaullEoOll HM) UV background. 
However, the spect ral energy distribution o f a starburst galaxy 
adopted from lBruzual A. & Chariot (1993, SB) leads to a rea- 
sonable model as well ( FechiieTetall2o"o* d). Therefore, another 
prediction is computed considering the results from the SB 
model for the z = 0.8643 system. The resulting metal line spec- 
trum is also shown in Fig. [3] The main differences in the pre- 
dicted lines are a decrease of the neon absorption (Ne m, Ne iv) 
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Fig. 2. The Keck data in the range z > 2.21. The data normalization is performed simultaneously with the estimation of the line 
parameters in course of the profile-fitting procedure. 



in case of the SB model, but slightly stronger Hei features. 
The predicted Oiv features are consistent with the observed 
data. However, the strong feature at 11 74 A is mainly due to 
Ov absorption from the z = 0.8643 system according to the 
HM model, while no O v absorption at all is predicted by the 
SB model. We find no alternative identification for this strong 
feature either proposed by the metal line system models or as 
part of the Lyman series of a strong, metal-free, low redshift 
system. Furthermore, no interstellar absorption line is expected 
at this wavelength. We conclude that also the energy distribu- 
tion of a starburst galaxy does not provide an optimum model 
of the z = 0.8643 system. Improvement of the photoionization 
model by adopting a different ionizing background is difficult 
since only few transitions are observed in the optical. The vast 
majority of lines constraining the model are loca ted in the UV, 
where the data suffers from low S/N (see also iFechner et alJ 
120061) . When we analyze the He n Lya forest taking into ac- 
count the metal line prediction, we will reconsider the lines of 
the z = 0.8643 systems with different strengths dependent on 
the adopted model (especially those of Ne iv and Ne v) and dis- 
cuss possible effects on the results. 

Absorption from galactic molecular hydrogen is expected 
to arise additionally in the FUSE spectral range. Templates for 
the H2 features with rest wavelengths in the far-UV includ- 
ing Lyman and Werner transitions and a detailed descripti on of 
how to identify and eliminate them can be found in McCandliss 
il2003l) . Since the proposed procedure requires the detection of 
a sufficient number of undisturbed galactic H2 features, it is 
inapplicable to our data. Instead, w e proceed as follows: The 
templates from M cCandlissI J2003I) are used to easily identify 
H2 lines, considering that they may be slightly shifted in ve- 
locity space. Then, we select unblended features in the FUSE 
spectrum and fit them with Doppler profiles to estimate the line 
parameters. The selected transitions and the fitted line param- 
eters are given in Table [2] All considered rotational transitions 
(ground state to J" — 0, 1 , 2, 3) can be modelled with a Doppler 
parameter b = 5.0 at a velocity v = -45.7 km s -1 . The de- 
rived column densities for every J" are also listed in Table |2] 



Table 2. Transitions of galactic H2 that are observed unblended 
in the FUSE spectrum. The atomic dat a (rest wavelength Ao an d 
oscillator strength /) are adopted from Ab g rail et alJ Jl993albl) . 
The column density is estimated assuming b = 5.0kms~' and 
v = -45.7kms _1 . 



transition 


-to (A) 


/ 


log A' 


L(3- 


0)R(0) 


1062.882 


0.0178 


16.00 + 0.87 


L(8- 


0)R(D 


1002.453 


0.0183 


17.83+0.25 


L(3- 


0)R(D 


1063.460 


0.0119 




W(0- 


0)R(2) 


1009.023 


0.0156 


14.72 + 0.30 


W(0- 


0)Q(2) 


1010.938 


0.0245 




W(0- 


-0)P(2) 


1012.170 


0.0055 




L(7- 


0)R(2) 


1014.977 


0.0190 




L(7- 


0)P(2) 


1016.458 


0.0102 




L(6- 


0)P(2) 


1028.106 


0.0105 




L(5- 


0)R(2) 


1038.686 


0.0165 




L(4- 


0)P(2) 


1053.283 


0.0090 




L(8- 


0)R(3) 


1006.413 


0.0158 


15.49 + 0.30 


W(0- 


0)R(3) 


1010.129 


0.0138 




L(6- 


0)R(3) 


1028.983 


0.0173 




L(5- 


0)R(3) 


1041.157 


0.0159 





Because of the few observed unblended lines and the poor data 
quality, the column densities, especially those of the lower rota- 
tional transitions, are very uncertain. Nevertheless, the derived 
values are roughly consistent wit h the column densities de- 
rived bv lRichter et aTlll200lll20oi for high- and intermediate- 
velocity clouds, respectively, towards extragalactic sources. In 
the final step, all transitions within the observed spectral range 
with J" = 0, 1,2,3 are included into our model. 

In the following we treat the H2 model like the predicted 
metal absorption lines. Unless explicitly noted we refer to the 
predicted metal line absorption as well as galactic H2 when 
talking about metal lines. Using either the profile-fitting proce- 
dure or the apparent optical depth method all additional lines 
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Fig. 3. Predicted metal line spectrum (smooth solid line) in comparison to the data (histogram-like solid line). The dotted line 
represents the error of the observational data. The predicted features are identified given the ion and the redshift of the system 
where they arise. 
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are considered in the analysis of the He n Lya forest. Besides 
the metal lines also the interstellar lines Cn A 1036 and Ari 
AA1048, 1067 are included. 

We roughly estimate how strongly the FUSE spectral range 
might be affected by lines of the higher Lyman series lines of 
the H i Lya forest at z < 0.2. The model spectrum derived in 
Fec hner et alJ (120061) is compared to the STIS UV data cov- 
ering the wavelength range 1230 - 1550 A used to constrain 
the models. Assuming all unidentified, strong features are low 
redshift Lya lines, the expected positions of the higher order 
Lyman series lines are computed. Most of the systems are not 
expected to contribute to the absorption in the FUSE spectral 
range since higher order Lyman series lines located redwards of 
the He n forest are detected only very weakly. However, three 
systems at z = 0.1240, 0.0809, and 0.0237 possibly affect the 
He n Lya forest. The Ly/3 feature of the system at z = 0.1240 is 
blended with Ly6 of the metal line system at z = 0.2140 both 
located at ~ 1 153 A. The corresponding Lyy and 6 features are 
expected at 1093.1 A and 1067.5 A, respectively. From the sys- 
tem at z = 0.0809 absorption is expected at 1108.7 A (Ly/3) 
and 1051.2 A (Lyy), while the Ly(3 features of the z = 0.0237 
systems should be located at 1050.0 A. As an be seen from 
Fig. |3 no remarkably high 77 values are detected at z ~ 2.456, 
2.460, 2.514, and 2.650, corresponding to the wavelength po- 
sitions discussed above. Therefore, biases in 77 due to higher 
order Lyman series lines of the low redshift H i forest should 
be negligible. 

The importance of metal line absorption in the FUSE spec- 
tral range is illustrated in Fig. |4] The upper panel shows the 
redshift distribution of 77 estimated by fitting Doppler profiles 
to the H i and He n Lya forest ignoring metal line absorption 
(for details of the analysis procedure see Sect. l5.1> . while metal 
line features have been considered in the lower panel. Lower 
limits indicate He n absorbers without detected H i counterpart. 
About 32 % of the apparent He n features have no H i counter- 
part when metal line absorption is ignored. The consideration 
of metal lines reduces the number of lower limits significantly, 
since several absorption features can be identified as metal 
lines. In this case only 25 % additional He n lines are needed. 
Data points that have been biased by metal line absorption, are 
marked in the lower panel of Fig. @] In total roughly 27 % of 
the measured values are affected. For 19 % of the data points, 
77 changes by more than A log 77 = 0.01 when extra absorption 
is taken into account. Considering only features of galactic H2, 
about 6 % of the data points would have been biased, i.e. the 
major contribution comes from the predicted metal line absorp- 
tion. There are no indications of unusually high column density 
ratios at the position of metal lines arising from the z = 0.8643 
system. 

The median column density ratio decreases from log 77 = 
1 .99 if metal lines are neglected to 1 .93 if they are taken into 
account. Consequently, the consideration of metal line absorp- 
tion has no severe influence on the general statistical properties 
of the 77 distribution. The column density ratio of a single ab- 
sorber, however, can be distorted up to an order of magnitude. 



5. Simulations 

Two different methods have been used so far to analyze the 
data o f the He 11 Lya forest towa rds HE 2347-4342. lKriss et alJ 
ll200ll) and lZhengeta"fl ll2004h) fit line profiles to the observed 
absorption features. The use of profile fits assumes the absorp- 
tion features to arise from discrete clouds in the IGM, which is 
certainly an oversimplification. The alternative method measur- 
ing the app arent optica l depth was also applied to the HE 2347- 
4342 data (Shull et al. 2004). In this case the column density 
ratio 77 is estimated by the ratio of the Hen and Hi optical 
depth per redshift or velocity bin. Since 77 is defined per bin, 
the assumption of discrete absorbers can be dropped in favour 
of a continuous medium, where the features are due to density 
fluctuations. A third approach fitting the optical Hi Lya for- 
est directly to the UV spectrum in order to derive 77 has been 
introduced by Fech ner & Reimersl J2006I) . 

In order to study possible effects of the analysis method on 
the results we create artificial datasets. By evaluating the arti- 
ficial spectra we investigate how the results are influenced by 
the applied method. For both methods we recognize problems 
which prevent the accurate reconstruction of the presumed 77 
value. In the following we briefly summarize the methodical 
procedures, then describe the generation of the artificial data, 
and discuss the implications of the profile-fitting procedure and 
the difficulties using the apparent optical depth method. 

5.1. Methods 

In preparation of the profile fit of the He 11 spectrum, we first 
identify all Hi Lya features in the optical data and estimate 
their parameters by fitting Doppler profiles. In the case of 
strongly saturated Lya lines higher orders of the Lyman se- 
ries up to Lye are considered as well to derive the param- 
eters. Then, the Hen Lya forest is fitted fixing the redshift 
and the Doppler parameter derived from the corresponding 
hydrogen lines. The latter constraint implies, th at we assume 
pure turbulent broadening. IZheng etalJ lE004b) made the at- 
tempt to derive the dominant broadening mechanism from the 
HE 2347-4342 data fitting unblended He 11 features with a free 
/^-parameter. They found a velocity ratio £ = £>He/£>H close to 
1 .0, suggesting turbulent broadening to be dominant. However, 
iFechner & Reimersl J2OO6) demonstrate that thermal broaden- 
ing is important for part of the absorbers and can lead system- 
atically to lower 77 values for high density Hi absorbers. 

The procedure of the apparent optical depth method applied 
to He 11 Lya forest data as described by Sh ull et alJ (J2004) will 
be summarized bri efly: The colum n density ratio is replaced by 
77 = 4 • THen/THi dlVIiralda-Escudel 1 993l more exactly the fac- 
tor 4 represents the ratio of the rest wavelengths /lo, h 1 / ^0, He n = 
1215.6701/303.7822), where the optical depths r = - \nF are 
measured per bin. In order to obtain physically reasonable val- 
ues, bins with normalized fluxes F above unity are masked out 
as well as bins with flux values below zero. Additionally, bins 
with flux values within 1 cr from unity or zero are not consid- 
ered, since they cannot be distinguished statistically from the 
continuum or zero, respectively. 
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Fig. 4. Comparison of the rj(z) distribution ignoring metal line and galactic H2 absorption (upper panel) and taking them into 
account (lower panel). Absorbers affected by extra absorption are marked as filled circles. Lower limits are based on Hen 
features without detected H 1 counterparts. Their Doppler parameters are chosen to be b — 27.0 km s . 



5.2. Creating artificial datasets 

Based on the statistical properties of the H 1 Lya forest as ob- 
served towards HS 1700+6416, a sample of lines is gener- 
ated. The column density distribut ion fu nction with (3 — 1.5 
is adopted from Kir kman & Tvtle Values in the range 

11.0 < logN < 18.0 are simulated. Our observed line sample 
yields j3 = 1.51 ± 0.05 for absorbers with 12.8 < logJV < 16.0. 
The Doppler parameter distribution is described by a truncated 
Gaussian 



dN(b) 
db 



A ■ exp 




( (b - bp) 
2o* 



2 \ 



b < b m \„ 



with bo = 27kms~ 1 , cr fe = 8.75 km s~', and b m [ a = lOkms -1 
as observed tow ards HS 1700+6416 in good agreement with 
iHu et al.l ( Il995l) . In addition, the parameters of the simulated 
line sample are correl ated by bmm - (1.3 • (\ogN - 12.5) + 
10.5) km s~' following Misawa (2002). The artificial spectrum 
covers the redshift range 2.292 < z < 2.555, which corresponds 



to wavelengths shorter than 1080 A in the FUSE spectrum. No 
number density evolution is considered along this redshift in- 
terval. The minimal distance between two lines is dz = 0.0001 
as observed towards HS 1700+6416. In the given redshift range 
the separation corresponds to ~ 20kpc. The resolution of 
R = 40000 and the signal-to-noise ratio of S/N ~ 100 are 
chosen to match the characteristics of high-resolution spectra 
taken with VLT/U VES or Keck/HIRES . 

The He 11 Lya forest is computed based on the artificial H 1 
lines using 77 = 80, the mean value found by Kriss et al .(2001). 
Considering a temperature of 10 4 K the Hen Doppler parame- 
ter is computed as 



?Hen = Jb 2 m -2kT 



\ m H inn, 



(1) 



The assumed temperature is coo ler than the average 2 • 10 4 K 
measured bv lRicotti et alJ J2000l) . but consistent with Z? m i n ,Hi = 
lOkms -1 . Test calculations with higher temperatures indicate 
no difference in the results. Since the IGM is expected to 
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Fig. 5. Distributions of the column density ratio 77 derived from 
the artificial data with a Doppler parameter consisting of a ther- 
mal and a turbulent part (upper panel) or a pure turbulent b- 
parameter (lower panel) applying a profile-fitting procedure. 
The dashed line indicates the presumed value 77 = 80. In the 
course of the analysis pure turbulent broadening is assumed in 
both cases. 



have a distinctive temperature distribution, a constant tempera- 
ture model is certainly an oversimplification, but serves to get 
an impression of the possible effects assuming pure turbulent 
broadening. For comparison we compute as well a He n dataset 
with £>Hen — bm- The resolution of R — 15 000 and the signal- 
to-noise ratio of S /N ~ 5 are chosen to match the typical values 
of real FUSE data. 

Since the artificial spectra are generated from statistical 
assumptions alone, they provide only a simple approach to 
investigate the applicability of the methods. More realistical 
spec tra based on hyd rodynamical simulations are presented 
bv lBolton et al.l (12006). They find that the column density ra- 
tio can be estimated confidently by a line profile-fitting pro- 
cedur e (but see discussion below and in iFechner & Reimerd 
2006). Furthermore, their simulations permit the investigation 
of different UV backgrounds, whereas the simple statistical ap- 
proach introduced here can provide only limited information 
about basic problems of the analysis procedures. 

5.3. Profile fitting analysis 

Having analyzed the artificial H i Lye forest by fitting Doppler 
profiles, the statistical properties of the sample are recovered. 
Under the simulated conditions our sample is complete down to 
log N ~ 12.0. Due to blending effects, the deduced Doppler pa- 
rameter distribution is broadened by 28 %. The He n Lyo- forest 
lines are treated the same way, using the derived H i parame- 
ters, with fixed line redshifts and ^-parameters. Features with 
Doppler parameters b < lOkms -1 are not considered. The re- 
sulting 77 values are shown in Fig. [5] Both the combined thermal 
and turbulent as well as the pure turbulent broadened sample 
show a scatter in 77 of about 0.5 dex. Generally, the line sample 
that contains thermal broadening has more low 77 values. We 
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Fig. 6. Recovered Hen column densities versus Hi column 
densities for combined thermal and turbulent broadening (up- 
per panel) and in the pure turbulent case (lower panel). The 
solid line represents the underlying correlation 77 = 80. For 
clarity, measurement errors are not plotted. 



find a statistical mean of log ?7thei- m = 1.77 ± 0.43 in compari- 
son to log7/t ur b = 1.86 + 0.37 in the pure turbulent case, both 
still consistent with the presumed value of 1.90309. The me- 
dian, which gives less weight to outliers, is log77 t herm = 1-86 
and log ?7tuib = 1 -90, respectively. 

There are various reasons for outliers lying outside the 
range log 77 = 1.903 + 0.500. Both models show only about 
3 % high values caused by line saturation and blending effects, 
which can, in principle, produce low as well high 77 values. In 
the sample with pure turbulent broadening we find 8 % low 
77 values while the thermal and turbulent broadened sample 
shows 14 % low values. Besides saturation and blending, the 
uncertainties of the parameters of weak Hi lines resulting in 
large error bars are a primary reason. In the case of combined 
thermal and turbulent broadening some He 11 lines are narrower 
than the adopted pure turbulent ^-parameter of H 1. This effect 
produces about 25 % of the low values. Therefore, if the widths 
of intergalactic absorbers are not completely dominated by tur- 
bulent mechanisms, low 77 values can be cau sed by the assump- 
tion o f pure turbulent broadening (see also Fec hner & Reimera 
2006). 

In order to find a strategy for avoiding artificial scatter 
in the 77 distribution, we directly examine the correlation be- 
tween Hi and Hen column density as presented in Fig. [6] 
The correlation fits well except for very low column densities 
(log A^h 1 < 12.0) where noise makes it difficult to derived accu- 
rate line parameters, and also for high column densities, where 
saturation of hydrogen lines plays a role (features are saturated 
when logA^Hi ^ 14.5 for b « 27kms _1 ). Since the redshift 
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and Doppler parameter of the helium lines are fixed during the 
fitting procedure, the column densities are derived correctly if 
the underlying assumptions are correct. In the case of combined 
thermal and turbulent broadening the He n column densities are 
significantly underestimated for log A^hi Z 13.0, caused by the 



assum ption of pure turbulent broadening. Fechner & Reim ers 
(2006) argued that the incorrect assumption of pure turbulent 
broadening affects only lines above this Hi column density. 
Thus, a stringent criterion to define a reasonable subsample 
would be 12.0 < logA^i < 13.0. Assuming pure turbulent 
broadening as a good approximation, 12.0 < logA^i < 14.5 
provides a less stringent constraint. 

II j -i 
Bolton et al. (2006) find a distribution similar to that shown 

in the lower panel of Fig. [6] Their lines samples contain only 

absorber with logA^i > 12.0 and show enhanced scatter for 

logA^Hi > 13.0 even in the case of a uniform UV background. 

Since they assumed turbulent broadening as well, the scatter 

most likely in dicates the effect of thermal line widths. However, 

iBolton et alJ ( |2006|) neither discuss the origin of the scatter nor 

the possible implications of the assumption of pure turbulent 

broadened lines. 

Applying the criteria derived above to the simulated 
datasets we find the statistical means log^erm = 1-85 + 0.21 
and log?7 tul b = 1.92 + 0.21 for the more stringent constraint 
(12.0 < logiVm < 13.0). In the pure turbulent case, the un- 
derlying j] value is recovered well, while it is still underesti- 
mated when thermal broadening plays a role. For both sam- 
ples the statistical 1 <x error has been reduced by a factor ~ 2 
in comparison to the total sample. Considering the sample 
12.0 < logA^Hi < 14.5 yields \og?] theim = 1.74 + 0.39 and 
log^tuib = 1.92 + 0.25, respectively. Here, the underestimation 
of j] for the data including thermal broadening gets worse and 
the scatter increases, while the mean value and 1 <x error for the 
pure turbulent data do not change significantly. 



5.4. Apparent optical depth method 

The same artificial datasets are analyzed applying the apparent 
optical depth method. The resulting distribution for a bin size of 
Av = 20 km s~ 1 , corresponding to 0.05 A in He n wavelengths, 
is shown in Fig.0 This method leads to a scatter in rj of about 
0.5 dex for both models as well. The high values are statisti- 
cally negligible, since they are only about 3 % for each sample. 
Again we find slightly more low rj values when thermal broad- 
ening is present (roughly 13 % in comparison to 10% in the 
case of pure turbulent broadening). This is plausible, since the 
He n optical depth in the wings of a thermally broadened, and 
therefore narrower feature is systematically lower than in the 
case of pure turbulence. In both models the underlying rj is un- 
derestimated. We find the mean values log 7/ t herm = 1-80 + 0.44 
and log?7 tur b = 1.84 + 0.41. With increased bin size the mean 
values get even lower (log 77 t h e rm = 1.75 ± 0.39 and log 77 tur b = 
1.77 + 0.38 for Av = 60kms~', corresponding to a bin size of 
0.20 A in Hen). Also the median underestimates the real 77 in 
the combined broadening case because of the excess of low val- 
ues, even though this effect is smaller. We find log J/therm = 1 -88 




Fig. 7. Distributions of the column density ratio r] derived from 
the artificial data with a Doppler parameter consisting of a 
thermal and a turbulent part (upper panel) or a pure turbu- 
lent ^-parameter (lower panel) using the apparent optical depth 
method. The dashed line indicates the presumed value 77 = 80. 
The bin size is 20kms~', corresponding to ~ 0.05 A in Hen. 
Bins with -2.0 < log Tm < —1.0 (see text) are denoted as filled 
circles, all others as open circles. The corresponding distribu- 
tions are indicated as histograms in the right panels. 

and log?/ tur b = 1.90 (bin size 20kms~ 1 ), in comparison to 
log?/therm = 1.81 and log77 tur b = 1.83 with 60kms~ 1 bin size. 

The correlation between low H i op tical depth, referred to 
as voids, and high 77 values as found by S hull et al.l 112004) can 
also be seen in the simulated data even though it is not present 
in the underlying spectra. For further investigation of this point 
it is instructive to regard 77 as the correlation between thi and 
THen- This is shown in Fig. [8] in the pure turbulent case for 
different signal-to-noise ratios of the He 11 data (100, 10, and 5, 
while S IN — 100 is constant for the H 1 spectra) and an example 
without any noise. 

At the given S/N of the Hi spectrum (= 100), optical 
depths down to logrni > -2.0 can be measured reliably, at 
least in principle. If only bins detached by more than 2 <x from 
unity are considered, the limit is logTm > -1.7. The scatter 
at low optical depths (e.g. for logTHi < —1.7 in the case of 
S IN = 100 for H 1 as well as He 11) is to due the noise. In order 
to avoid a contamination by noise-affected bins, the constraints 
defining the subsample should be chosen carefully. As can be 
seen from Fig. [8] the detection limit of the Hen optical depth 
depends on the signal-to-noise, of course. In case of S /N — 5 
the restriction is log Tjien > -0.65, or log rnen > -0.3 for the 
2 <r constraint. 

Even the analysis of the spectra without noise leads to scat- 
ter around the underlying relation. Furthermore, bins with an 
optical depth above log thi > —1.0 (indicated by the vertical, 
dotted line in Fig. [8} yield Hen optical depths clearly beneath 
the expected values. This flattening is due to line saturation, 
since saturated absorption features become broader and more 
smeared out at finite resolution. The difference in the profile 
does not only affect the saturated core but as well the wings 
of the line. Consequently, the flattening starts at relatively low 
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Fig. 8. The Hen optical depth versus Hi optical depth for 
simulated data with pure turbulent broadening. The specific 
signal-to-noise of the artificial helium dataset is given in each 
panel. The signal-to-noise ratio of the hydrogen spectrum is 
S IN — 100. The upper panel is based on artificial datasets 
with no noise at all. The resolution of the simulated spectra 
are /?hi = 40000 and /?Hen = 15 000, respectively. Each point 
represents a bin of 20kms~ 1 . The solid lines indicate 77 = 80. 



hydrogen optical depths. Thus, a direct comparison of the ap- 
par ent optical depths above a c ertain limit is not possible. 

ISavage & Sembachl Jl99ll) give a detailed description of 
the method and applicability using the apparent optical depth. 
They emphasize the possibility of finding line saturation in 
the analysis of line doublets. Hidden saturation can be seen if 
the apparent column density N a (v) oc T a (v)/f Aq of a doublet 
pair differs with respect to the considered component. Since 
line saturation is important in the whole FUSE spectrum, the 
apparent optical depth method is not applicable to bins with 



thi > 0.1 if the underlying 77 is roughly 80. Adopting a 1 cr 
detection limit of thi > 0.01 there is 1 dex of hydrogen opti- 
cal depth left, allowing the apparent optical depth to be applied 
reasonably. 

IFox et all d2005l) investigated the effect of noise on the ap- 
parent column density. They found that the apparent optical 
depth method will likely overestimate the true column density 
when applied to data with low S /N. Since the FUSE spectra 
are very noisy (S /N ~ 5), the results obtained by an apparent 
optical depth method have to be considered carefully. 

The bins selected by -2.0 < logrei < —1.0 are marked 
in Fig.Qand the overlay in the right panels shows their 77 dis- 
tribution. Both the samples (pure turbulent and combined line 
widths) contain virtually no low 77 values and the fraction of 
high values is below 5 % in each case. As expected, the mean 
values approach the underlying 77 and the scatter decreases. 
We find log^therm = 1-88 ± 0.24 and log77, urb = 1.90 + 0.24, 
respectively, with a bin size of 20kms _1 . Fitting a Gaussian 
A ■ exp(-(log 77 - log ?7o) 2 /2 a 2 ) to the log 77 distribution of the 
subsample leads to log 770 = 1 .90 + 0.01 and cr = 0.237 + 0.009 
in the case of pure turbulence. In the case of combined ther- 
mal and turbulent broadening we find log 770 = 1.89 + 0.16 and 
cr = 0.213 + 0.012. In comparison to the mean value of the 
sample without any noise, which is 1.89 + 0.07 in the pure tur- 
bulent case, these values demonstrate that the scatter in the 77 
distribution is mainly due to the noise level of the data. With 
the present FUSE data (S/N ~ 5) an 77 distribution with a scat- 
ter of roughly 13 % can be expected following the simulations. 

5.5. Comparing profile fitting and apparent optical 
depth 

The investigation of artificial spectra by fitting Doppler pro- 
files or applying the apparent optical depth method exhibits dif- 
ferences of the robustness and applicability of the diagnostics. 
Concerning the apparent optical depth, line saturation in com- 
bination with the detection limit given by the noise leads to a 
relatively small interval of 1 dex in the hydrogen optical depth 
(0.01 < thi < 0.1), where the method can be applied reliably. 
Besides limited resolution noise is the most crucial limitation, 
since measuring the apparent optical depth per velocity bin is 
directly sensitive to the noise level, while the impact of the as- 
sumption of pure turbulent line width is less significant (the 
deviation is roughly 1 %). 

Since position and line width are fixed, noise has less influ- 
ence using the standard profile-fitting method. Therefore, even 
saturation is less problematic. Hydrogen starts to be saturated 
at column densities logA^Hi ^ 14.5 and normally higher order 
Lyman series lines can be used to estimate the H 1 line param- 
eters accurately in case of strong saturation. Thus, even sat- 
urated Hen features may be evaluated reasonably. However, 
in the crowded Lycr forest the number of components may be 
ambiguous leading to larger uncertainties. Another point is, in 
analyzing the simulated data, we know that the correct profile 
function is applied since the datasets are created using it. When 
dealing with real data, the observed features can be described 
only approximately. 
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Table 3. The rj values recovered from the simulated data us- 
ing the profile fitting or apparent optical depth method, respec- 
tively. The input value is 80, i.e. log rj = 1.90309. The first 
columns show the values for the profile-fitting analysis con- 
straining the sample to 12.0 < logA^i < 13.0 (selected) and 
without any restriction (unselected). The columns at the right 
hand side give the results obtained with the apparent optical 
depth method (binning 20km s~ ) with -2.0 < logrni < —1.0 
(selected) and for the total sample (unselected). 



profile fitting apparent optical depth 
(log;;) median (log 77} median 



turbulent 


selected 


1.92 


±0.21 


1.92 


1.90 


±0.24 


1.89 




unselected 


1.86 


±0.37 


1.90 


1.84 


±0.41 


1.90 


thermal 


selected 


1.85 


±0.21 


1.87 


1.88 


±0.24 


1.88 




unselected 


1.77 


±0.43 


1.89 


1.80 


±0.44 


1.88 



A sample restricted to lines that can be clearly identified 
and are insensitive to effects of thermal broadening is expected 
to yield the most reliable results. The selected lines have to be 
strong enough to be fitted unambiguously but are below the 
sensitivity limit for thermal line widths (12.0 < logA^i < 14.5 
or even less), since a line sample under the simplistic assump- 
tion of pure turbulent broadening will underestimate the true 77 
values. According to iFechner & Reimersl J2OO6) features with 
log A^hi 5 13.0 remain unaffected by effects of thermal broad- 
ening. Therefore a more stringent selection would be 12.0 < 
logjVm < 13.0. 

Table[3]lists the recovered 77 values and their statistical 1 cr 
errors using either the profile fitting or the apparent optical 
depth method. The median reproduces the underlying value 
of log 77 = 1 .903 within an accuracy of more than 93 % in 
all cases. The statistical mean underestimates the column den- 
sity ratio fitting line profiles, particularly if therm al broadening 
is pres ent. This is consistent with the results of iBolton et alJ 
(2006), who come to the conclusion that the median gets closer 
to the underlying value, although they do not discuss possible 
reasons for this finding. Furthermore, the scatter is reduced by 
more than 40 % if only the selected subsamples are considered. 

Following these results, we expect to find a scatter in log 77 
of roughly 20 % by analyzing the observed spectra without any 
further selections. Applying the apparent optical depth method 
to a subsample with 0.01 < thi < 0.1 leads to a reduced scat- 
ter (13 %). For the profile-fitting procedure the scatter can be 
reduced as well down to about 1 1 % if only absorbers with 
12.0 < logA^Hi < 13.0 are selected. This stringent constraint 
minimizes the effect of systematic underestimation of 77 if ther- 
mal line broadening is present. 

6. Results and discussion 

Having an idea of the advantages and limitations of the analysis 
methods, we present applications to the observed data and dis- 
cuss their implications. In contrast to the simulated data, the ob- 
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Fig. 9. Redshift distribution of the column density ratio 77 for 
absorbers with 12.0 < logA^i < 13.0. The dashed line indi- 
cates the median value log 77 = 1.89. 

served spectra contain additional absorption features of metal 
lines and galactic molecular hydrogen absorption bands. The 
modelling of this additional absorption, as described in Sect. 
|4j is considered in the following analyses. In the course of 
the profile-fitting analysis, the metal line parameters are fixed 
during the fitting procedure. Using the apparent optical depth 
method, bins with a metal line or H2 optical depth r met > 0.05 
are omitted. 

6.1. Profile fitting analysis 

In the He n redshift range covered by FUSE (z > 2.29) we find 
326 H 1 Lya lines in the Keck spectrum. Due to the FUSE de- 
tection gap and terrestrial airglow lines, the number of H 1 lines 
that are supposed to be detectable in He 11, reduces to about 300. 
In addition to the He n absorbers, which are identified in H 1, we 
have to add several He n lines without any detected H 1 coun- 
terparts. These absorbers are about 25 % of all fitted He n lines. 
Furthermore, we include the metal absorption lines described 
in Sect. 0] The resulting distribution of the column density ra- 
tio 77 with redshift is shown in the lower panel of Fig.|4] 

Because of the high fraction of lower limits, i.e. Hen fea- 
tures without H 1 counterpart, the simple statistical mean of the 
absorbers with detected H 1 and He n (log 7/ = 1 .54 + 1 .01, me- 
dian 1.72) is expected to be biased. Usin g the Kaplan-Meier es- 
timator (e.g. Feigelson & Nelson 1985) to derive a mean value 
including the lower limits yields log77KM = 2.13 + 0.08. The 
median of the total sample is log 77 = 1.93 in good agree- 
ment with the results from the studies of the Hen Lyor for- 
est towards HE 234 7-4342 dKriss et al.l200lUShull et al.l2004 
IZheng et all2004bl) . revealing 77 * 80. 

The values given above are estimated considering the total 
line sample. As we argue in Sect.|5] more confident results are 
expected to be found if the line sample is restricted. Therefore, 
we compute the mean column density ratio considering lines 
with 12.0 < logA^Hi < 14.5 and yield log 77 = 1.63 + 0.82 
(median 1.72). Since for all Hen features without Hi counter- 
part the lower limit of the hydrogen column density is less than 
10 12 cm 2 , no limits are included in the selected line sample. 
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Applying the more stringent selection criterion (12.0 < 
log Mi i < 13.0) results in log 77 = 1.85 + 0.70 (median 1.89). 
The redshift distribution of rj for this subsample is presented 
in Fig. [5] Even for the smallest sample the scatter of log rj 
is roughly 40%. This is significantly more than expected by 
the simulations. A possible explanation would be the redshift 
evolution of the column density ratio. Due to the evolution of 
the sources and absorbers that compose the properties of the 
UV background, p is expected to change with redshift (e.g. 
iFardal etal .1 ^998). Our artificial spectra are generated free 
from any evolution. Thus, the real data might mimic a higher 
overall scatter due to redshift evolution. We estimate the col- 
umn density ratio in different redshift bins of the most strin- 
gently selected sample and find indeed an increase of the mean 
value with redshift. For z < 2.50 we yield log// =s 1.7, while 
log 77 ~ 2.0 is estimated at z > 2.60. However, the scatter in 
each redshift bin is about 20 - 40 %, i.e. part of the scatter in 
the total redshift range might be due to evolution, but a signif- 
icant fraction cannot be explained neither by statistical scatter 
nor by redshift evolution. A Kolmogorov-Smirnov test between 
the results from the artificial data and the observations yields a 
probability of less than 1 % for each redshift bin that both sam- 
ples are based on the same distribution function. This leads to 
the conclusion that the UV background is indeed fluctuating, 
even though the results from the profile-fit analysis are insensi- 
tive to the scales of the 77 variation (an attem pt to derive these 
scales is made bv lFechner & Reimersl20 06). 

For investigating the behaviour of 77 with the hydrogen col- 
umn density, we concentrate on absorbers with log Mi 1 < 15.0. 
The H 1 features stronger than this threshold are associated with 
the LLS in the observed redshift range. Since LLS are believed 
to arise from material in the outer parts of galaxies, their col- 
umn density ratio does not probe the ionizing conditions in 
the IGM but in a galaxy itself. These absorbers are excluded 
from the sample for the following investigation. Furthermore, 
we consider only absorbers detected in both H 1 and He 11 with 
reasonable column density uncertainties cr(logMii,Hen) < 1-0. 
The distribution of log/7 with log Mi 1 is shown in the upper 
panel of Fig. ^3 A clear trend is seen, that log 77 decreases with 
increasing H 1 column density. We compute the Spearman rank- 
order correlation coefficient yielding r s = -0.64. A linear fit to 
the data log 77 — a ■ log Mi 1 + b leads to a — -0.62 + 0.05 and 
b = 9.67 + 0.64. 

The analysis of the artificial data (Sect. [5} reveals that 
a correlation arises if a sample of lines broadened by ther- 
mal and turbulent processes is analyzed under the assumption 
of pure turbulent broadening. The effect is especially notice- 
able for the column density ratio of absorbers stronger than 
logiV(Hi) > 13, as illus trated in Fig.|6|and d iscussed in the pre- 
vious section (see also lFechner & ReimerslEo06l) . In this case 
we find a Spearman rank coefficient of about -0.25, while no 
correlation is found if the assumption of pure tur bulent broad- 
ening is correct (r s « 0). Indeed, the results of IZheng et alJ 
(2004b) imply t hat the fine width is dominated by turbulent 
broadening, but iFechner & Reimersl J2006I) have shown that 
thermal broadening cannot be neglected completely. 

Applying the more stringent constraint (12.0 < logMii < 
13.0, presented in the lower panel of Fig. fTot to the observed 
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Fig. 10. Distribution of the column density ratio 77 with hydro- 
gen column density Mii (both logarithmic). Only absorbers de- 
tected in Hi and Hen with cr(log Mi 1, Hen) < 1-0 are shown. 
The upper panel presents the total sample with log Mi 1 < 15.0. 
The lower panel is a zoom-in to 12.0 < log Mi 1 ^ 13.0. 
The dashed line represents the median value for the more con- 
strained line sample log 77 = 1.89. 

data yields a correlation coefficient r s = -0.33 for the observed 
sample, which is only half of the value for the larger sample. 
The linear fit shows now a flatter slope a = -0.54 + 0.13 and 
b = 8.75 + 1.66. Also in this case, the Spearman rank coeffi- 
cients indicate a slight correlation for the artificial sample with 
fractional thermal broadening (r s =s -0.2) and no correlation 
for the pure turbulent broadened sample. Since we derive a 
stronger anti-correlation from the observed data compared to 
the artificial spectra including thermal line widths, the conclu- 
sion that 77 is slightly larger in voids might be correct. 

6.2. Apparent optical depth method 

Considering the limitations discussed in Sect. [5] we apply the 
apparent optical depth method to the observed data. In a ddi- 
tion, w e use the metal line system models of IFechner et alJ 
(2006) with the modifications discussed in Sect.|4]to omit bins 
which are affected by metal line absorption in the FUSE or the 
Keck data. Fig. II II shows the resulting 77 distribution for both 
the whole data and the subsample selected by -2.0 < log r Hl < 
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Fig. 11. Distribution of 77 derived form the observed data of HS 1700+6416 using the apparent optical depth method with a bin 
size of Av = 20kms~ 1 omitting bins affected by metal line absorption. The selected subsample (-2.0 < logT H i < —1.0) is 
indicated as circles and its 77 distribution is the narrower histogram in the right panel. The dashed line in the left panel represents 
the median value of the subsample log 77 = 1 .90. 



-1.0. The subsample consisting of 501 points (54 % of the total 
sample) only contains values in the range of 1.0 < log 77 < 3.0 
with a mean value of (log 77) = 1.90 ± 0.30 (median 1.90). A 
consideration of the total sample yields log 77 = 1 .65 + 0.55 and 
a median of 1 .7 1 . 

Fitting a Gaussian to the log 77 distribution of the subsam- 
ple yields log 770 = 1.903 + 0.013 and cr = 0.287 + 0.010 
corresponding to FWHM = 0.676 + 0.023. The width of the 
distribution is increased by more than 30 % in comparison to 
the artificial data. Thus, we may conclude, that in addition to a 
scatter in 77 due to noise and saturation effects a physically real 
variation of 77 is present. 

Fi gures ITU and ITU show that 77 is decreasing with redshift. 
The presented 77 values are based on the effective optical depth 
of H 1 and He n, respectively, per Az = 0.05 bin, where censored 
pixels have been neglected. In the resulting evolution 77 is rising 
from log 77 = 1.70+0.32 at z ~ 2.33 up to 1.98+0.32 at z ~ 2.68, 
consistent with the results from the profile-fit analysis. This 
behaviour of the He 11/H 1 ratio is expected, since it traces the 
evolu tion of the ionizing sources. According to iFardal et"a"fl 
(1998), the observed evolution may be mainly based on a pop- 
ulation of quasars with spectral indices a 1.8. 

An important source of error regarding the redshift evolu- 
tion of 77 are the uncertainties in the placement of the continuum 
in the FUSE spectrum, since the He 11 optical depth and thus the 
value of 77 directly depends on the position of the continuum 
level. This is also true for the estimate of the column density 
in case of line-profile fits. As discussed in Sect.[3]the choice of 
the extinction curve, which is highly uncertain in the EUV, has 
a strong impact on the continuum level. Especially, if the con- 
tinuum is extrapolated to a too low level at short wavelengths 
(e.g. the dotted curve in Fig.[JJ, the apparent optical depth and 
consequently the derived 77 gets too small. Since the extinction 
law is wavelength-dependent, the uncertainties increase with 
decreasing wavelength. Additionally, the inaccurate definition 
of the zero flux level at short wavelengths results in systemati- 
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Fig. 12. Evolution of 77 with redshift. The dotted curve repre- 
sents the expected values if the ionizing sources are quasars 
with a spectral index of a = -1.8 as modelled bv lFardal e t al. 
<199§. 



cally underestimated He n optical depths leading to low 77 val- 
ues. Thus, it might be an artifact from the continuum definition 
and/or zero flux level estimation that the observed evolution 
of log 77 appears to be steeper than predicted by th e model of 
IFardal et alJ lfl"998). However. Eheng et alJ d2004bh find a very 
similar evolution towards HE 2347-4342, a line of sight with a 
lower extinction (E(B - V) = 0.014). 

A correlation between the H 1 optical depth and the He n/H 1 
ratio is still present in the subsample. The Spearman rank-order 
coefficient for this correlation is r s = -0.53. However, also 
the subsamples of the artificial datasets show slight correlations 
leading to r s ~ -0.36 for both broadening mechanisms. The 
slope describing the observed sample (a = -0.66 + 0.04) is 
roughly twice as steep as those from the simulated datasets with 
uncertainties of the same order of magnitude. This is consistent 
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Fig. 13. Distribution of the spectral indices a (used as f v oc y a ) 
for the observed sample (derived from the appa rent op tical 
depth result; thin line) and the whole sample from lTelfer et alJ 
(2002, larger bin size, thick line). The observed spectral in- 
dices are derived from the measured 77 values using the correla- 
tion illustrated as Model A2 in Fig. 10 of lFardal et ail dl998l) . 
Different bin sizes are used for a better presentation. 



with the results from the profile-fit analysis discussed in detail 
in the previous Section. We therefore confirm the conclusion 
that part of the correlation between H 1 voids and high 77 values 
might be real. 

6.3. The ionizing background 

Summarizing the results obtained by applying the apparent op- 
tical depth method and the profile-fitting procedure, we find a 
mean value of the column density ratio log 77 = 1 .90 (corre- 
sponding to 77 « 80). The observed scatter is larger than ex- 
pected from the analysis of artificial data. Part of the enhanced 
scatter is possibly due to evolution of the column density ra- 
tio with redshift, which has not been incorporated in the arti- 
ficial data. Measuring 7/ and its scatter in smaller redshift bins 
suggests as well that part of the scatter must be attributed to 
physical effects beyond statistical noise. 

The origin of the scatter must be related to fluctuations in 
the He 11 ionizing radiation since the photoionization rate of 
H 1 can be assumed to be uniform (a mo re detailed argumen- 
tation is given by e.g. Bolt on et al.ll2006l) . Spatial variation at 
the He 11 ionization edge may be due to the wide r ange of spec- 
tral indices of the energy distribution of quasars dZheng et alJ 
1 19971 ITelfer et al.ll2002[ IScott et alJEool . IShull et alJ d2004 
and lZheng et alJd2004bl) compared the estimated column den- 
sity ratios towards HE 2347-4342 to the distributi on of the 
the sp ectral indices a of 79 quasars presented by ITelfer et ail 
(2002). Both found an excess of softer radiation. Fig.ll3lshows 
the distribution of a derived from the observed subsample ap- 
plying the apparent optical depth method. The 77 values are 
converted into spectral indices using the relation between the 
He II /Hi ra t io and the source's spectral index presented by 
iFardal et all dl998l their Fig. 10, Model A2). As a mean spec- 



tral index, we find (a) _ —1.99 ± 0.34. In comparison to the 
total sample of ITelfer et alJ d2002l) . we confirm the excess of 
softer spectral indices. One should keep in mind that the con- 
version from 77 into a is based on a phenomenological model, 
which certainly contains additio nal sources of erro r. 

Following the arguments of IShull et alJ J2004), it is possi- 
ble to convert the 77 values directly into effective spectral in- 
dices. The effective spectral index a e s describes the radiation 
an absorber is directly exposed to rather than the radiation once 
emitted by the sources, i.e. it denotes the spectral index of the 
filtered radiation of the sources. Assuming a temperature of 
T = 10 4 3 K and the equality of the local spectral indices at 
1 a nd 4 Ryd, respectiv ely, leads to 77 = 1 .70 • 4~ ffclr (for details 
see IShulletalJ (|2004) and references therein). Applying this 
relation the a distribution is broader and shifted to even softer 
indices: {a e g) = -2.78 + 0.50. The results from the profile- 
fitti ng method lead to s imilar numbers. 

iBolton et"aTld2006l) have studied the observed spatial vari- 
ation of the column density ratio arising from a fluctuating 
UV background. According to their model, the fluctuations 
are based on variations in the Hen photoionization rate due 
to properties of the emitting quasars, such as the spatial dis- 
tribution, lumin osity function, and the spectral energy distri- 
bution (see also lFardal etaDll998l) . lBolton et al.1 J2006I) find 77 
valu es scattered oyer 2 de x roughly consistent with the results 
from lZheng et alJ d2004bl) . but failed to reproduce a significant 
number of points with 77 < 10. As discussed in the previous 
Sections, the results from the line profile-fitting procedure are 
likely to be biased due to thermal line broaden ing when con- 
sideri ng high de nsity absorbers (see also iFechner & Reim ers 
2006). Although Bolto n et al.l d2006l) do not address this prob- 
lem their Fig. 4 reveals an increased number of lines with a low 
column ratio for increasing Hi column density. Considering 
only absorbers in a restricted column density range eliminates 
most of the points with very low 77 (see Figs.l9land ll It . We con- 
clude that there are no absorbers exposed to radiation harder 
than those emitted by QSOs. Thus, the column density ratios 
inferred from the data of HS 1 700+6416 support the model 
proposed bv lBolton et alJ d2006l) to explain the spatial fluctua- 
tion of the column density ratio. 

Apart from scatter, the average value of the column den- 
sity ratio can be roughly related to the dominating sources 
of ionizing radiation. The inferred value (77 « 80) is in 
very good agreement with the results towards HE 2347 -4342 
dKriss et alJboOlt IShull et alJEool IZheng et alJl2004bl). It is 
consistent with the UV background of Haardt & Madau (2001) 
based on QSOs only. However, there is growing evidence that 
there are contributions to the ionizing background from other 
sources than quasar s ("e.g. [Heap et al.|2p00||Smette et alJ2002t 
lAguirre et allho04l IBolton et alJl2006t iKirkman et alJl2005l) . 
The excess of soft spectral indices found from Fig.ll3lsuggests 
that part of the absorbers are indeed exposed to softer radiation. 

Our data reduction may still be in need of improvement 
concerning the zero flux level (see Sect. This is in par- 
ticular true for shorter wavelengths. In order to exclude pos- 
sible biases due to the reduction problems, we estimate the 
mean column density ratio for wavelengths longer than 1086 A 
(z > 2.58), where the data reduction appeal's to be correct. We 
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Fig. 14. Evolution of the He n opacity with redshift. In addi- 
tion to the opacities measured in this work (filled circles), val- 
ues measured towards five additional QSOs with detected He n 
absorption are sho wn. The solid curv e represents the relation 
Teen <* (1 + z) 35 JFardal et alj fl998) fitted to the data points 
Z < 2.75. The dashed line indicates an overall fit to the data 
resulting in Teen oc (1 + z) 62±04 . 



find a median of log 77 = 2.01 adopting the profile-fitting re- 
sults and a slightly lower value (x 1.97) using the apparent 
optical depth method. These values (77 « 93 . . . 102) are con- 
sistent with predictions for t he UV bac kground including also 
young star forming galaxies (lHaardt & Madaul200ll) and in the 
same range as predicted by the models of iBolton et alJ J2OO6) 
at the same redsh ift based on a modified version of the UV 
background of Mada u etaD (Qj 



6.4. Redshift evolution of the He 11 optical depth 

Fig.ll4lshows a compil ation of the measured He 11 opacity Thp tj 
from the literature (Davidsen et al. 1996; Anderson et al. 1999; 
iHean etalJEoffll IZheng et alJl2004alht iReimers et alJl2005a^ 
together with our measurements from the FUSE spectrum of 
HS 1700+6416 (filled circles). Estimating the mean normal- 
ized flux per redshift bin, pixels affected by metal line ab- 
sorption are excluded. We find a mild evolution from THen = 
0.83 + 0.24 at z = 2.30 - 2.35 to 1.06 + 0.18 at z = 2.70 - 2.75 
with a slight dip (0.74 + 0.34) at z ~ 2.45. These values are con- 
sistent with the opacity mea surements towards HE 2347-4342 
in t he same redshift range ftZheng et all2 004b). 
iDavidsen etal.llfl99i 

estimated THe n = 

1.00+0.07 at <z> = 

2.4 using a low resolution spectrum taken with the Hopkins 
Ultraviolet Telescope (HUT). As can be seen from Fig.ll4lour 
values are systematically lower by about 25 %. Tests with the 
data uncorrected for metal line absorption suggest that this dif- 
ference is partly due to the additional absorption by metal lines 
raising the opacity by roughly 10%. Additionally, differences 
in the continuum leve l may account for the h igher effective op- 
tical depth derived bv lDavidsen et alJ ([1996). However, the val- 
ues are consistent on the 1 <x level. 

The s olid line in Fi g. ITU represents the relation THen K 
(z + l) 3 ' 5 (Far dal et alJl9 98) providing a suitable fit to the data 



at z < 2.75. The value of the power law exponent y + 1 is cho- 
sen according to the results from H 1 Lya forest statistics (e.g. 
IPress et alJll993t iKim et alJl2002l) . Deviations from this rela- 
tion at higher redshift are commonly interpreted as the evidence 
for the tail end of the epoch of He 11 reionization which com- 
pletes at z ~ 2.8 ( e.g. IReimers et afll 19971 iKriss et alJEoOll 
Zhen g et ali20 04b). At z > 3 very high values of the He 11 opti- 
cal depth can be detected. A simple overall fit to all data points 
would lead to T He n = 4.0 • 1 0~ 4 (z + 1 ) 6 ' 2±a4 (dashed line in Fig. 
114-1 . This fit should not be interpreted as a reasonable model for 
the evolution of THen- 

7. Summary and conclusions 

We have presented far-UV data of the quasar HS 1700+6416 
taken with FUSE, which is the second line of sight permitting 
us to resolve the He n Lya forest. The data are of comparable 
quality (S /N ~ 5, R ~ 20000) to those of HE 2347-4342 and 
cover the redshift range 2.29 < z < 2.75. In this redshift range, 
no strong variations of the He n opacity are detected. The evo- 
lution of the effective optical depth is consistent with (z + l) 3 ' 5 , 
i.e. this line of sight probes the post-reionization phase of He n. 

The column density ratio 77 has been derived using line pro- 
file fits and the apparent optical depth method, respectively. A 
preliminary study of simple artificial spectra created on the ba- 
sis of the statistical properties of the Lya forest including or 
neglecting thermal line broadening, respectively, reveals the 
shortcomings of the standard analysis methods. The reasons 
are noise, line saturation, and effects due to the assumption of 
pure turbulence if thermal broadening contributes to the line 
widths. We find that the profile-fitting procedure leads to re- 
liable results only for absorbers with Hi column densities in 
the range 12.0 < logA^i < 13.0. The apparent optical depth 
method is only valid in the range 0.01 < Tm < 0.1. Otherwise, 
77 would be underestimated in the case of strong H 1 absorption. 
Furthermore, a scatter in 77 of 10 - 15 % is expected even if the 
underlying value is constant. 

In order to avoid systematic biases due to additional absorp- 
tion in the He 11 Lya forest, a model of the met al line feat ures 
expec ted to arise in the FUSE spectral range jFechner et alJ 
2006! has been included in the analysis. We have also con- 
sidered features of galactic H2 absorption. 27 % of all fitted 
Hen lines are affected by metal lines, in case of 19% of the 
lines the derived 77 value changes significantly, i.e. by more than 
A log 77 = 0.01. Additionally, the required number of Hen ab- 
sorbers without detectable H 1 counterpart is reduced by 30 % 
if metal lines are taken into account, and the average 77 value 
slightly decreases. Although the consideration of metal line ab- 
sorption does not distort strongly the statistical properties of 
the resulting 77 distribution, individual values could have been 
biased by up to an order of magnitude. 

Photoionization models of the absorbers showing associ- 
ated metal lines provide an independent estimate of 77. For these 
systems the He 11 column density can be computed based on the 
photoionization model for the metal line features. Comparing 
the 77 values directly measured from the profile-fitting proce- 
dure and inferred from the metal line system modelling leads in 
principle to consistent results. A more detailed discussion will 
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be given in a future paper dFechner & Reimersll2006l) which 
will also discuss the absence of the He n proximity effect to- 
wards this QSO. 

For the redshift range 2.58 < z < 2.75, where the spectrum 
appears to be free of artifacts due to the reduction process, the 
data reveals log rj « 2.0 on average. This value (77 » 100) indi- 
cates a contribution of galaxies to the UV background at these 
redshifts, consistent with current results from studies of the H 1 
opacity and the metallicity of the IGM. 

The scatter of r\ is larger than expected compared to our 
analysis of artificial datasets. Therefore, we infer that the UV 
background might really fluctuate, even though our present re- 
sults are insensitive to amplitude and scale of these variations. 
The main limiting factor for a quantitative estimate of the fluc- 
tuations is the high noise level of the data. Converting the val- 
ues of the column density ratio into spectral indic es, we con- 
firm t he appar ent excess of soft s ources as found bv lShull et alJ 
(2004) and lZhengetai1 (l2004b). Our results suggest, that the 
apparent correlation between the 77 value and the strength of 
the Hi absorption may be an artifact of the analysis method 
and does not reflect physical reality. 

Since we have shown that the assumption of pure turbulent 
line broadening will lead to systematical errors, the next step in 
the analysis of the resolved He 11 Lya forest should be to avoid 
this assumption. The apparent optical depth method applied to 
an appropriately restricted sample should lead to unbiased re- 
sults independent of the dominant line broadening mechanism. 
A bin by bin comparison as performed by the apparent optical 
depth method provides little potential for taking into account 
a potentially dominant thermal broadening, while the profile- 
fitting procedure could drop the assumption of pure turbulence. 
However, obtaining reasonable fits would be more challenging 
due to the low S /N of the data. This would provide an indepen- 
dent strategy to estimate the IGM temperature. 
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